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Introduction. 

Quantum evaporation of black holes exhibits fascinating connection between gravity, 
quantum mechanics and thermodynamics. In order to understand this phenomena better, 
and certainly in order to move on towards solution of bigger questions, like the back- 
reaction, the issue of the final state, or the possible evolution from pure to mixed state, 
(otherwise known as the "loss of quantum coherence"), it would be usefull to examine 
different paradigmas and different techniques that have been used in the past for study of 
black hole evaporation. Historically the first has been Hawking's use of mode mixing and 
Bogoliubov coefficients [1] . Close to this is a technique utilizing properties of the squeezed 
states in black hole spacetimes, which received some attention recently [2]. Next, there 
is a puzzling and not yet fully understood discovery of a relationship between the flux 
of the outgoing radiation and the trace anomaly [3], which has been used most recently 
in studies of dilaton models [4]. Another approach explored in seventies has been in the 
physical analogy with the Klein paradox [5]. But most people consider as the most superior 
the method of thermal Green functions, introduced by Gibbons and Perry [6]. This is just 
the method that we will explore here, but our starting point will be a comprehensive study 
of the analytic properties of Green functions, done by Fulling and Ruijsenaars [7]. We 
will pragmatically follow their path, by computing first a number of relevant quantities for 
scalar field in flat spacetime at zero and at the finite temperature. Next we move on to 
curved spacetime and discuss thermal properties of Rindler and Schwarzschild spacetime. 
This methods allows for easy and clean evaluation of the temperature of the black hole as 
seen by the stationary observer outside the horizon. 
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Scalar field in fiat spacetime 



The basic relations for a massive scalar field are 



= / d 4 x 
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u(k) = +\J k 2 + m 2 . 



= ; [$(t,x),7r(t,$)]=i5(x-$) . 



Suppose now that we want to compute the two point corellation function 



(*(a;)*(y)) • 



The two interesting choices for the quantum state in which this expectation value is to 
computed are vacuum and the thermal state. We will consider them now in turn. 



Vacuum expectation values 



If we choose to compute expectation values in the vacuum state, it follows that, 



<0|$(x),$(y)|0)= / dk t 



_ ik-(x-y)-ik°(x°-y ) 



2u{k) 



Using Cauchy formula, 

i f ... m 



f( z o) = w~ I dz 

Zm JC(zn) 

we can rewrite this as, 



C(z ) z- Zq 



(o|»W«te)|0> = (-0/g I (t o^ t)a ■ < 2 > 

The correspodence between these two expressions is not unique, as the later integral 
is multivalued. The integrand has two poles, at k° = ±u>(k), and the final result depends 
on how does the contour of integration go around the poles. This choice of contour is 
equivalent to the choice of boundary condition for the corellation function. The residua at 
the two poles are (— 1) / (2k^ ole ) , and the two basic contributions are 

' e"*"**)* , and —L-e™™* . 



2u{k) 2w(fc) 
So, the general result will be, 



(0\*(x)Mv)\0)c = I j^y 3 e^" ^ H)^ 



. .UO ' 

where C stands for the contour of integration. The 2-point correlation function in Eq. 
(2) corresponds therefore to one particular choice for the contour of integration in (2), or 
to one particular choice of the boundary conditions for the solutions of Eq.'s (3-4) below. 
Different choice of contour leads to different vacuum expectation value in Eq. (1). 
Different vacuum expectation values have been historically defined as follows: 

G + (x,y) = 9(x° — y°){0\<&(x)<&(y)\0) (positive frequency Wightman function); 



G (x,y) = 8(y° — x°)(0\3>(y)$>(x)\0) (negative frequency Wightman function); 
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iG(x,y) = (0| [<&(#), <E»(y)] 1 0) (Pauli - Jordan function); 



G W (x,y) = (0\{$(x),$(y)}\0) (Hadamard function); 



Gr(x, y) = (—)9(t x — t y ) G(x,y) (retarded Green function); 



G^(x, y) = 6(t y — t x ) G(x, y) (advanced Green function); 



iG F (x,y) = (0\T($(x)$(y))\0) (Feynman propagator); 



All of these 2-point functions are usually called Green functions although four of them 
actually obey the homogeneous equation. There is a number of useful and easy to prove 
relations that are obeyed by these functions: 

Problem 1. Show that, 



iG = G + - G~ ; = G + + G~ ; iG F = G+ + G~ . 



when g e {G+,G-,iG,GM}. 

Problem 3. (a) Show first that, d? t (6G) = 5G + 6G; then, (b) show that the other 
three vacuum expectation values are true Green functions: 



Problem 2.. Show that, 



(-d 2 x + m 2 ) Q(x,y) = 



(3) 



(-d 2 x + m 2 ) g(x,y) = ±5(x-y) 



(4) 
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when Q G {G_r, Ga, Gp}. Plus sign applies for the first two, minus sign for Gp. 

Further claim is that solution to any of these equations has the following integral 
representation: 



For Q G G~ , iG, G^}, the right hand side of this equation should vanish; when 

Q G {Gr, Gai Gp}: the integrand of the Fourier transform on the right hand side should 
be ±1. Both can be accomplished using the multivaluedness of the complex integral (5), 
as the following excersizes will show: 

Problem 4- Draw all seven contours of integration for defined Green functions, and 
discuss boundary conditions. 

Problem 5. Since first four functions obey the same equation, they are distinguished 
by different boundary conditions, and similarly for the other three: 

(a) Determine boundary conditions at ±oo for four solutions to the homogeneous 
equation, and justify names and notation for G^. 

(b) Determine boundary conditions for the three solutions to the inhomogeneous equa- 
tion (4) by adding a source term J(x)$(x) to the original Lagrangean. 

Let us now fix x y ^ x cts ct parameter, and consider all these Green functions 
as functions of time. It is immediately obvious that, by virtue of its defining expression, 
G + (t) may be analytically extended only to the upper half-plane, and G~ (t) only to the 
lower. Further straigthforward examination uncovers the following beautifull structure on 
the complex plane z = (t, r): 




(5) 



which obeys, 




6 



Problem 6. Show that there is an analytic function Q(z) which coincides with the G + 
on the upper half-plane, and with G~ on the lower half-plane. Along the t axis Q has 
a branch cut for t > \x\. The jump accross the cut is given by the Pauli- Jordan Green 
function: Q{t + ie) - Q(t - ie) = G(t), for e > 0. 

Problem 7. The projection of G to the imaginary axis is called the Euclidean, or the 
Schwinger Green function: Q(z) — > t _».o Ge(t). Show that: 

(a) Ge obeys the Euclidean version of the inhomogeneous field equation for a Green 
function, 

(d 2 TT + V 2 + m 2 )G E (r,x) = -5(r)S ( - 3 \x) ; 

(b) Ge and Gf are related by the so-called Wick rotation: t — > it. 

Results of these excersize may be stated as follows. There is Green function Q{z), 
which in upper half-plane represents the analytic continuation of G + (t), and in lower half- 
plane it is the analytic continuation of G~(t). Along the real axis Q has a branch cut for 
t > \x\. The jump accross the branch cut is given by the Pauli- Jordan Green function, that 
is, by the canonical commutator for the scalar field <fi. On the imaginary axis Q reduces 
to the Euclidean, or Schwinger Green function Ge(t). Unlike all the others mentioned so 
far, Ge is a true Green function: it obeys the Euclidean version of the field equation with 
the point source. Ge(t) by itself has its own analytic continuation to the real axis: its is 
the Feynman propagator, Gp(t), which obeys the inhomogeneous field equation with the 
point source. This continuation is usually called the Wick rotation. 

This description, first given so eloquently by Fulling and Ruijsenaars, is the culmina- 
tion of our understanding of the Green functions. 

We can now compute explicitly these functions. The most efficient strategy is to com- 
pute G ± first, by evaluating the integral above with the appropriate boundary conditions. 
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Other Green functions may be then computed from G ± , by adding or subtracting the two 
residua according to their contours, and/or the algebraic relations between the functions. 
Results of the excersizes below should be compared to the calculations in textbook of 
Bogoliubov and Shirkov [8]. 

Problem 8. Compute G + for a massless scalar field in flat spacetime. Since there is 
no free parameter in the Lagrangean there is no characteristic correlation scale and this 
Green function describes a long range correlations. 

Problem 9. Compute G + for a massive scalar field in a flat spacetime. This Green 
function does have a characteristic correlation length, given by the finite mass of the 
particle. In the perturbation theory this would lead to the Yukawa potential. 

Problem 10. Now compute all the other Green functions for a massless scalar field in 
flat spacetime. 

Problem 11. What can we conclude from these Green functions about the physics of 
a massless scalar field? 

Problem 12. Compute all the Green functions for a massive scalar field in flat space- 
time. 

Problem 13. What can we conclude from these Green functions about the physics of 
a massive scalar field? 

Thermal expectation values 

Instead of the vacuum, we now consider the thermal state, 

p = "£\E)e-< 3E (E\ . 

E 

The 2-point correlation function is, 

(<S>(xMy)) = J2e- PE (EM*My)\E) . 

E 
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It is straightfoward to show that the later can be written in the path-integral form. For 
this we should imagine that (3 is the total interval in the imaginary time within which the 
"evolution "takes place: (3 = Tf — T{. Separating this interval into the infinitesimal ones, 
and inserting the unit operators represented through complete sums of the coordinate or 
momentum eigenstates will lead to the formal expression, 



Here C($) stands for the class of the fields over which the path integral is taken. In this 
case they are those with the periodic boundary condition in the imaginary time: 



This kind of reasoning is widely used, but it usually leaves some ambiguity as to what 
is exactly the relationship between the fields and the Physics in the real and imaginary 
times. Practice has shown what is the safe way to interpret this relation in many cases, 
but the situation often gets ambiguous when the curved spacetime is introduced. It is in 
this methodological respect that the construction of Fulling and Ruijsenaars appears to us 
so effective. We will now repeat the preceding analysis from the vacuum case to determine 
the analytic structure of the Green functions in the case of a free field in the flat spacetime, 
but at the finite temperature T = f3~ l . 

The two Whightman functions are defined as follows: 




C($) = {$(r,r) | ®{r,T i + P) = Q(r,T i )-,T i <T<Tf = T i + P} . 



G+(t, x, y) = Z- 1 Tr{e"^$(t, x)&(0, y)} ; 



n 



Gp(t, x, y) = Z- 1 Tr{e-^$t (0 , f)} 



n 



The usual interpretation is that is an amplitude to create particle at (0, y) and anni- 
hilate it at all in the a thermal bath (rather than vacuum). 
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Let iV+ denote the number of particles with positive energy u n > 0, and N n number 
of antiparticles with the same energy. The partition function is given as, 

Z= J2 e _/J £« (JV » +JV ' ;> " = JKl-e-^™)" 2 

{N+,N~} 

The power of two appears here if we are dealing with a complex, non-hermitean field. If 
$ = <E>t, then, 

z = H(i-e- /3 ^y 1 . 

n 

As in the vacuum case, where (0|oa|0) = (0|a+at|0) = 0, we have Tr{paa} = TrjpaV} 
0. However, instead of (0|a^a|0) = 0, we have the finite occupation numbers for thermal 
state: 



so, that, 



(a n c4) =5 nrn (l-u^ ^ log z\j =S nm (l-e 1 , 



(ala n ) = -1 + (a n al) = (e^« - l) 1 



Using the expansion for the quantum field with time dependent piece explicitly written 
out (the spatial piece depends on the choice of the spatial coordinates and topology) , 

fn(x 



we find the following expressions for the Whightman functions: 



fn(x)f*(y) 



2U 7 



-iuj n t 



+ 



I — e -/3ui n e /3w„ _ I 



e ~iuj n t e iu n t 



+ 



Let us now look at the analytic coninuation of these functions. The first step is to 
establish their relationship to the Pauli- Jordan function: 
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Problem 14- Check that still, 



Since exponents with both signs are present in the expressions for G^, if we now 
substitute z = t + ir for t, we can see that neither of the function is analytic on the whole 
half-plane, either upper or lower. However, observe that if we fix r > 0, for u n — > oo, 

G+ ~ e w " T , Gp ~ e w " (T -^ ) . 

Similarly, for r < 0, the leading pieces are, 

Therefore, Gp may be extended from the real axis to the strip r G [— /?, 0], and Gp 
to the strip r G [0, /?]. 
Furthermore: 

Problem 15. Observe that, 

G+(r-(3)=Gp(r) ; Gp (r + (3) = G+ (r) . 

Using this result it is easy to see that one can do analytic continuation to the whole 
plane, using strips of the width (3 in imaginary direction as Weierstrass circles. One 
therefore arrives on the following analytic structure on the complex plane: there is a 
complex function Gp(z,x), periodic in imaginary direction with the period (3. Within 
each strip [r, r + Q has two representations, one through Gp, another through Gp, 
but with identical values. This realization within one strip is repeated within all other 
strips. Along the real axis, for t > \x\, Gp has a branch cut. The jump accross the cut 
is given by the canonical commutator, that is, by the Pauli- Jordan Green function at 
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the zero temperature. Due to the periodicity in r this branch cut is therefore present at 
r = 0, ±/3, ±2/3, ±3/3, .... 

Let us now consider the restriction to the imaginary axis, the so-called thermal Eu- 
clidean function, or thermal Schwinger function, 



Gp(z, x) ^t^o Gp (r, x) = 9(x) 



+ 



/3oJ n _ I I _ e -U> n (3 



In zero temperature case Ge was an even function, Ge(—t) = Ge(t). Here we find 
the same result: 

Problem 16. Show that G|(-r) = Gf (t + /?) = G|(r). 

Problem 17. Show that those two properties, periodicity in r and even nature of Gp, 
combined together make this function periodic with a period (3/21 

As people say, Gp is reflected around r = n/3/2. This is then repeated through all 
the other ,3-wide stripes. 

Problem 18. Draw the analytic structure of thermal Green function Gp(z) in the 
complex plane z = t + ir. 



Rindler spacetime 

Let us now consider the case of a flat spacetime as seen by an uniformly accelerated 
observer. The trajectory of such an observer moving in X direction is given by, 

t = A' 1 sinh AT , X = A' 1 cosh AT . (6) 

Here A e [0, oo] is constant proper acceleration, and T e (— oo, oo) is his proper time. For 
future convenience let us introduce the dimensionless affine parameter A = AT, and the 
"radius" a = A~ l . 
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If we consider the set of all the observers with all possible values of constant proper 
acceleration A we see that pair (A, a) may be used as a new set of coordinates in place of 
(£, X). The line element may be written as, 

ds 2 = -a 2 d\ 2 + da 2 + dl\ , 

where the last term stays for the length in two unchanged coordinates. As it is well known, 
due to limitations on the speed to which observer may accelerate to, these coordinates cover 
just the two wedges in the (i, X) plane, within the straight lines X = ±t. But within these 
wedges the spacetime looks flat in new coordinates (A, a). 

Consider now the thermal expectation value as measured by such an observer: ( A, a) . 
Without any calculation, we know from the preceding discussion the analytic structure of 
this function, in particular we know that it must be periodic in imaginary proper time 
S : T -> S with a period /3. 

Let us compare this with the vacuum expectation value measured by the inertial 
observer using Minkowski coordinates (t,X). This would allow us to compare the two 
physical pictures as seen by the two observers. The punch line is, of course, that those two 
pictures correspond to the same physical situation, as we shall see shortly. To facilitate 
such comparaison we only need to do a very innocent step, just to perform the change of 
coordinates in one function to those naturally used for another. We must also take care of 
the proper dimensions, as it will become apparent. 

We will start with the vacuum expectation value as calculated by the Minkowski 
observer. Given two points, (ii, Xi, Yi, Z\) and fe, X%, Y2, Z%), the vacuum expectation 
values must depend only on As 2 = —At 2 + AX 2 + A/ 2 , due to the Lorentz invariance. 

After substitution of the Rindler coordinates, viz. Eq. (6), we obtain function that 
depends only on a finite line element expressed as, a 2 + a\ — 2a\(i2 cosh(Ai — A2) + A/|- 
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One can now make a series of simple observations. 

First, note that upon change to the purely imaginary Rindler time, T — > iS, equiv- 
alently, A — > zu, the line element becomes, ds 2 = da 2 + a 2 da 2 . This looks like the line 
element in polar coordinates, and there will be no singularity as a if a is indeed a 
polar angle. But in that case it must be periodic, with a period 2tt. The corresponding 
imaginary proper time S must then also be periodic, with a period 2%/A. That is, we use 
A = AT generalizes to A + ia = A(T + iS), so that a = AS. In that case every function of 
S must also have the same period, or fraction of it. In particular, we can see that for the 
vacuum expectation value evaluated above the Grand vacuum Green function function G 
has the form, 

Q{a\ + a\ — 2aia 2 cos(ai — a 2 ) + dl\) , 

along the imaginary axis. This is the same periodicity as the thermal Green function in 
flat spacetime at finite temperature A/(2tt)\ 

This type of the reasoning is used very often: periodicity of the metric in imaginary 
time is interpreted as the presence of thermal features. However, to make reliable statement 
one has to investigate the analytic structure of Green functions. 

The elementary observation is that Q(t,X), being the vacuum expectation value, is 
holomorphic within the strip At < AX, i.e., for As 2 > 0. But this line element is invariant 
under the change of coordinates, hence, this expectation value, as seen by the observer 
using Rindler coordinates (A, o) has the same holomorphic strip and the same periodicity 
in imaginary direction as what he would call the thermal expectation value. 

As a final confidence test, one can that the following statement holds: 

Problem 19. Show by a direct coordinate substitution from (6) that the equation 
for the Euclidean Green function in Minkowski space becomes the equation for Euclidean 
Green function in Rindler coordinates. 
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To summarize: After transformation to Rindler coordinates, <3(t+ir, Ax) has the same 
periodicity, and the same analyticity strip as Gp(T+iS, Ax) for (3 = 2% /A; Euclidean Green 
functions in Minkowski and Rindler coordinates are solutions of the two equations related 
by the simple change of coordinates; the same must be said about the boundary conditions 
they obey; thus, it is the same Green function; analytic continuation of Euclidean Green 
function from one streep to another along the imaginary axis uniquelly defines Q over the 
whole complex plane; therefore, 

g(t + iT,Ax)=g 27T/A (T + iS,Ax), (7) 

in the entire complex plane. 

On basis of that statement we conclude that while observer using Minkowski coordi- 
nates thinks all the time that he is computing averages with respect to the vacuum state, 
the Minkowski vacuum, observer using Rindler coordinates, the accelerated observer, no- 
tices that expectation values for a system (the scalar field) in the same state look like 
thermal expectation values, computed for a finite temperature field theory at tempera- 
ture A/(2tt). In other words, Rindler observer sees Minkowski vacuum as thermal density 
matrix, temperature being proportional to its proper acceleration. 

This is the main result that we wish to demonstrate here. Let us just mention that an 
apparent puzzle of seeing pure state as thermal may be explained by invoking horizons for 
the Rindler observer. Similarly, there is no problem of a finite state. With an unlimited 
supply of energy, Rindler observer may indefinitely maintain its constant acceleration. 
Once it stops accelerating, horizons disappear, and the two observers agree that they are 
seeing the vacuum state, as both are inertial observers. 
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Schwarzschild black hole 



This case offers considerably more puzzle, but thanks to its similarity with the Rindler 
spacetime it is easy at least to take a first step and establish the existence of thermal 
features. 

The standard form of the line element is, 

d8 2 = -(l-^)dt*+ ^ . + r 2 dQ 2 . (8) 
V r J 1 -r s /r 1 v ' 

The Schwarzschild radius is rs = 2GM, M being the mass of the black hole. 

Consider now an observer keeping himself at a fixed radius r above the horizon of the 
black hole. In order to do so he must accelerate away from the hole, with an acceleration 
just equal to the gravitational acceleration provided by the hole. Thus, such an observer 
is Rindler observer. 

The line element (8) may easily be transformed into the Rindler form in case of an 
observer not too far from the horizon. Let us have r = rs + Sr = r^(l + x), with x > 0. 
The line element becomes 

ds 2 = -xdt 2 +rl— . 

x 

This can be transformed to the Rindler form —a 2 d\ 2 + da 2 if we define a = 2r,g£ 1 / 2 , and 
consequently, 

A = . 

2r s 

Upon extension to the complex plane we should find, as before, the periodicity of the 
Green functions in imaginary dimensionless time a, with period of 2n. The corresponding 
period in the imaginary component of the Schwarschild time is then (3 = Anrs- Thus, the 
observer at a constant distance above the horizon of the Schwarzschild black hole sees the 
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thermal bath at temperature, 



What is this temperature of? If the observer that maintains constant distance from 
black hole is Rindler observer, then inertial observer is obviosly the one in a free fall 
towards the hole. Indeed, by equivalence principle, its vacuum state must be identical to 
the vacuum of the observer at infinity, which does not detect gravitationaly the existence 
of the black hole. Thus, if initial state of the quantum field is vacuum with respect to 
the empty space, once black hole moves in, (or we move towards the hole), we perceive 
quantum fluctuations in that state as decohered thermal fluctuations. 

And why do we call this temperature the black hole temperature? One can show, 
just as in the case of a Rindler observer, that horizon may be understood as the source of 
the black-body radiation associated with this temperature. By overall energy conservation 
this radiation must lead to the decrease in the mass of the hole, hence to its evaporation. 



Conclusion 

As said before, there are several ways to exhibit the thermal features of black holes. 
There is probably no the best way, and there is no need for one. The full understanding of 
the black hole evaporation most likely will require shifting between different paradigmas 
and methods of calculation. But the method reviewed here is the one that is most firmly 
rooted in the quantum field theory. It clarifies what does it mean to have finite temperature, 
and allows for a simple and reliable calculations. Thus, the use of thermal Green functions 
ought to be usefull in future studies of the final state of the evaporation, and the questions 
about the "loss of the quantum coherence." The main obstacle is how to proceed with 
what is essentially a non-equilibrium configuration, and how to take into account the 
back-reaction. 
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